Gradient Catastrophe and Fermi Edge Resonances in Fermi Gases 
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Any smooth spatial disturbance of a degenerate Fermi gas inevitably becomes sharp. This phe- 
nomenon, called the gradient catastrophe, causes the breakdown of a Fermi sea to multi-connected 
components characterized by multiple Fermi points. We argue that the gradient catastrophe can be 
probed through a Fermi edge singularity measurement. In the regime of the gradient catastrophe 
the Fermi edge singularity problem becomes a non-equilibrium and non- stationary phenomenon. 
We show that the gradient catastrophe transforms the single-peaked Fermi edge singularity of the 
tunneling (or absorption) spectrum to a sequence of multiple asymmetric singular resonances. An 
extension of the bosonic representation of the electronic operator to non-equilibrium states captures 
the singular behavior of the resonances. 



1. Introduction The FES (Fermi edge singularity [T]- 
[3]), observed as a power law peak in the absorption spec- 
trum of X-rays in metals more than 70 years ago, is one of 
the most prominent and well understood quantum many- 
body phenomena caused solely by Fermi statistics. 

FES also has been demonstrated in tunneling experi- 
ments [4]-[7]: a sudden switch-on of a contact potential 
due to a change in the capacity of the contact in tun- 
neling causes a power law dependence of the tunneling 
current on the bias voltage: I(V) ~ y-^a+ka Here 
6 = ttcl is the scattering phase of the ensuing potential 
and k is the number of scattering channels. In the case 
of an attractive potential (a > 0) the current peaks at 
the Fermi edge. 

The physics of the FES is explained by the phe- 
nomenon of the Orthogonality Catastrophe [9 : the state 
of a Fermi gas (Q f \ after a localized potential is suddenly 
switched on, is almost orthogonal to a state of the un- 
perturbed Fermi gas Their overlap vanishes with the 
level spacing A as a power law (fi'|fi) ~ (A/EF) ka . 

The FES acquires new features in the non-stationary 
regime due to the gradient catastrophe. A gradient catas- 
trophe is a hydrodynamic instability observed in many 
classical (and, recently, in atomic) systems. In Ref. [10] 
it has been shown that a quantum analog of the gradient 
catastrophe also takes place in a degenerate Fermi gas. 
The ballistic propagation of macroscopical packets and 
fronts in Fermi gases inevitably enters the gradient catas- 
trophe regime, where the initially smooth fronts develop 
large gradients and undergo a shock wave phenomenon: 
packets overturn as shown in Fig [T] . The observation 
of non-stationary phenomena in Fermi gases may not be 
easy since electronic times are too short, but does not 
seem impossible. From the theoretical viewpoint a non- 
stationary FES reveals important (and new) aspects of 
the Orthogonality Catastrophe. Both catastrophes are 
caused solely by Fermi statistics and therefore their in- 
teraction is of interest. 

The physics of non-stationary processes in Fermi gases 
is in its infancy. Ref. [10] discusses oscillatory corrections 
to the Orthogonality Catastrophe in a non-stationary 
regime: the overlap of the state of a shaken-up Fermi 
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FIG. 1: Time progression of the Fermi edge. The front over- 
hangs, giving rise to three Fermi (see inset) edges (Pf3 > 
Pf2 > Pfi) between the trailing and leading (x+) edges. 



gas with a propagating packet (both before and after 
the shock). We also mention Refs. [11] which cast non- 
stationary Fermi gases in the context of integrable non- 
linear waves. 

In this paper we study the FES in a non-stationary 
regime before and after the quantum shocks had oc- 
curred. We show that each shock introduces two ad- 
ditional Fermi edges, each edge causes an additional res- 
onance peak schematically depicted in Fig. [2] In general 
settings we predict that a smooth wave packet passing 
through a tunneling contact shows a sequence of singu- 
lar pulses in the tunneling current. Here we discuss only 
one shock, though the extension to multiple shocks is 
straightforward. 

Non-stationary FES in its generality consists of sev- 
eral different regimes some of them are difficult to study. 
In this paper we consider a regime away from turning 
points. There the Fermi surface changes much slower 
than the Fermi scale. In this case the problem becomes 
quasi-stationary and can be solved in two steps. First 
we neglect the time dependence of the Fermi points lead- 
ing to FES problem with a non-equilibrium population 
of electronic states. Once this has been done, we may 
treat the Fermi-edges as slow time dependent variables. 
We argue that they obey Riemann equation. 

The stationary non-equilibrium FES problem has been 
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extensively studied in Refs. jT2] and more recently in 
Refs. [T3UT5] , so we could have borrowed the results. We 
suggest instead a simple compact approach to the non- 
equilibrium FES which captures the singular parts of the 
resonances. This approach extends the representation of 
the electronic operator by means of Bose fields to non- 
equilibrium states. In a related paper [15] we formulate 
the non-equilibrium FES problem as a matrix Riemann- 
Hilbert problem with an integrable kernel. 

2. Tunneling in a non- stationary regime. We con- 
sider the following situation: 

(i) A Fermi gas is in contact with a localized resonant 
level (a quantum dot). It is initially uncharged and pro- 
vides no scattering to electrons. When an electron tun- 
nels to the dot, it suddenly charges the dot, switching-on 
a small potential H — >> H' = H + U localized at the dot 
[8]. We assume the potential to be weak, such that the 
scattering phase \S\ < tt/2 and therefore \a\ < 1/2. As- 
sume no further interaction, no dissipation, ignore spin 
and channels. 

(ii) A semiclassical electronic front or a packet - a state 
with a spatially inhomogeneous density matrix g has been 
initially created in a Fermi gas as it is shown in Fig. [I] 

This setting can be realized in different ways. For ex- 
ample, a smooth potential well (namely, of spatial extent 
much larger than the Fermi length) centered away from 
the dot is applied to the Fermi gas. A large number of 
electrons is trapped in the well. Then the well is suddenly 
removed. The electronic packet propagates towards the 
dot, eventually hits it, facilitating tunneling. We may 
consider two conductors with different chemical poten- 
tials suddenly brought into contact. The conductor with 
smaller chemical potential is in contact with a dot. A 
front of electronic density propagates toward the smaller 
chemical potential of the dot. 

Alternatively one may apply a time dependent voltage 
through a point like contact separated from the dot [16] . 

A state created in such processes is a Fermi coherent 
state 
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obtained by a unitary rotation of the ground state |0) of a 
Fermi gas. A function Pq(x) <C Pf is an initial momenta 
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FIG. 2: A schematic plot of the tunneling current for a > 
(left panel) and a < (right panel), solid lines show computed 
power law asymptotes. Dashed lines interpolate between res- 
onances, a < displays a single peak after the shock, while 
a > displays two peaks. 



of the packet. The field ip is a chiral canonical Bose field 
related to the chiral part of electronic density 



e ikx Pk, 



pk 



= ^2 c i c P+k, (1) 



where p and p + k are electronic momenta close to the 
Fermi momentum Pp. We assume that Pq(x) is smooth 
on the Fermi length scale H\VPq\/Po <C Pf- This condi- 
tion justifies the semiclassical analysis described below. 

The tunneling current is given by the golden rule [TJ [8] . 
In units of a tunneling amplitude I(oj)\noj=eV-\-E F reads 

/(cj,t)cxRe / e lUJT G(t+-,t--)dr, (2) 
Jo 2 2 

G(t u t 2 ) = (n\e iHt2 ce iH, ^- t2) c^ e - iHtl \Q) (3) 

Here c = ^ k e~ zkx °Ck and xo is the position of the dot. 

Since all the physics is concentrated at the Fermi edge 
a knowledge of the dispersion at the edge is sufficient 



e p = E F + v F (p- Pf) 



(p-Pf) 2 
2m 



(4) 



In the literature the parabolic part of the dispersion is 
routinely ignored. In this approximation our effect dis- 
appears [T7j . 

We evaluate G(ti 1 t 2 ) in the regime where the typical 
time of tunneling r = t 2 — ti > is much smaller than 
the time it takes for the packet to change. This approx- 
imation does not allow us to compute the broadening of 
the resonance at the frequency range 7 ~ vp(y Po/K) 1 / 2 , 
but it captures the power law shoulders of resonances and 
their dependence on time t. 

Under this assumption during the short time of tunnel- 
ing the energy dependence of the Fermi velocity and scat- 
tering phase 5 caused by the potential U can be dropped 
in some interval |e — Ep\ <C A at the Fermi edge, where 
the cut-off A (typically of the Fermi scale) is assumed to 
be larger than vpPo an d h/r. This amounts a shift of 
energy levels after scattering downwards by a constant 
amount a (in units of level spacing): e p — > e p — a. 

In Ref. [18] it has been shown that the vertex op- 
erator e aif implements a shift of momenta such that a 
perturbed Hamiltonian and perturbed states are H' = 
e -a<p(x ) He a<p(x ) and |^ = e M*o)|f^ Then Green's 

function reads 

G(t u t 2 ) = (^|c(t 2 )e- a ^^ ' £2) e a ^ ( ^' £l) c t (ti)|^). (5) 

This formula is standard. The only difference is that the 
density matrix does not commute with the Hamiltonian 
and therefore the process is not-stationary - the Green 
function and a current depend on t = 1/2 (ti + 1 2 ). 

3. Gradient Catastrophe: Riemann equation for Fermi 
Gases. We demonstrate the gradient catastrophe on the 
evolution of the Wigner function - a simpler object than 
The Wigner function describes occupation in phase 
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space: 



n F (x,p, t) 



J ("let 



(x+ -,t)c(x 



,t)\n) t 



4 dy 



(6) 



We assume that the front is plane or radial, such that 
the dynamics is essentially one dimensional and chiral. 

Semiclassically, the Wigner function is equal to 1 in a 
bounded domain p < P F (x,t) of the phase space (p,x) 
- the Fermi sea - and vanishes outside the Fermi sea 
n F (x,p,t) « Q(P F (x,i) — p). This form is valid as 
long as the gradients of the spatial dependence of the 
Fermi momentum P F (x,t) are small [20]. The shape of 
the initial Fermi surface is given by the density matrix 
P F (x, 0) = P F + P${x). The support of the Wigner func- 
tion is the area below the Fermi surface in Fig. [I] 

How does the Fermi surface change in time? It does 
not, if one neglects dispersion of the Fermi gas, i.e., treats 
the velocity v p = de p /dp as a constant: the front trans- 
lates with the Fermi velocity without changing its shape. 
It does change, in a dramatic fashion, if the dispersion in 
Q (no matter how small) is taken into account. 

The Wigner function (for a dispersion e p = p 2 /2m) 
obeys the equation 

(d t + v p V)n F (x,p,t) = 0, v p =p/m (7) 

The solution of this equation 

n F (x,p,t) = n F (x — v p t,p, 0) w 6 \Po(x — —t) — p) 

\ m / 

(8) 

shows that a moving Fermi momentum Pp(x,t) obeys a 
hodograph equation 



P F (x, t) = P (x - P F (x, t) jm ■ t) 



(9) 



This is Riemann's solution of Euler's equations for hydro- 
dynamics of a compressible one-dimensional fluid, also 
called Riemann (or Riemann-Burgers, or Riemann-Hopf) 
equation 

d t P F + VE F = 0, E F (x, t) = P F (x, t)/2m (10) 

The Riemann equation leads to shock waves: the veloc- 
ity of a point with momentum P F (x) is P F (x)/m: higher 
parts of the front move faster. The front gets steeper, 
and eventually attains an infinite gradient - a shock 
at some finite time. After this moment the Riemann 
equation has at least three real solutions, P F s(x,t) > 
P F 2(x,i) > P F ±(x,i), confined between two turning 
points x-(t),x+(t), the trailing and leading edges respec- 
tively (Fig{l]). 

This phenomena is the gradient catastrophe. Any 
smooth disturbance of the Fermi surface eventually ar- 
rives to a point where the Fermi surface acquires infinite 
gradients, and then becomes multi- valued between mov- 
ing turning points x±. We focus on that region, namely 
the region where all electronic states with energy below 
Efi and between Ep2 and E F s are occupied. The Fermi 
distribution acquires at least three, or more edges. 



5. Slowly evolving Fermi edges Away from turning 
points we can employ the Whitham averaging method 
known in the theory of non-linear waves [19]. The 
method has been applied to electronic systems in [TO] . 
It is based on a separation of scales between the slowly 
varying Fermi points and fast oscillations of the elec- 
tronic states. In short, the Whitham method suggests 
treating the slowly changing Fermi edges as constants 
while computing Green's function (J3|, and then to in- 
clude the motion of the Fermi edges in the final result. 
Motion of the Fermi edges is determined by the Riemann 
equation (TToj) . This approach is valid away from turning 
points |2l| ). It can be justified mathematically using an 
integrable non-linear equation for Green's function ob- 
tained in [11]. We omit the mathematical justification of 
this procedure since the approximation of slowly moving 
Fermi-edges is physically obvious. 

In the shock region X- < xo < x + we must evaluate 
the current ([3| over a state with three Fermi edges, where 
electrons occupy states below E F i and between E F 2 and 
E F s (see Figj3|, and then treat the time dependent edges 
as solutions of the hodograph equation Q . This is a non- 
equilibrium (and stationary FES problem) addressed in 
Refs. [IMS]. 
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FIG. 3: Unoccupied electronic states are between edges E F i 
and E F 2 and above E F 3- 



6. Fermi edge resonances Now we are ready to 
present our result: the current at a frequency close to 
the edges <C £i (hu — E Fi ) <C v f Pq : in units of cut- 
offs reads: 



t) ex \ fajj — E Fi 



3a —2eia 



3& 



-2£na 
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nm 



(11) 



where ^ = (-1) <+1 , i = 1, 2, 3 and = E Fi -E Fj . The 
Fermi energies all depend on space- time according to the 
hodograph equation We suppress this dependence in 
the notations. 

The last factor in this equation gives an overlap be- 
tween states before and after the shake-up - an Orthog- 
onality Catastrophe formula for multiple edges 



(12) 



where A is a level spacing. It is a power of Cauchy de- 
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constructed 



out of lower and upper edges of occupied bands. These 
formulas extend the result in [THS] to multiple edges (see 

also P2H33]). 

If the potential is attractive a > (a common case) 
the current features a peak at E F i (almost zero bias) 
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and an additional resonance at E F3 with a power law to 
the right of the edges. Current is suppressed at the edge 
Ep2- If potential is repulsive a < a peak may appear 
at the edge Ep2 with a power law to the left to the edge 
Figg 

Apart from additional resonances, the unique features 
of the shock region is the value of the exponent and the 
time-dependent amplitudes. Outside of the shock region 
where there is only one Fermi edge or at a larger energy 
\Epi — Epj \ <C \uo — Ep\ < A where the fine structure of 
Fermi edges becomes negligible, the current is given by 



E F 



where 



the standard formula [TJ [2] I(uj) oc (uj 
only the edge Ep depends on time. 

7. Bosonic representation for non- equilibrium states 
Techniques developed in Refs. [T2HT5] allow computing 
the details of the FES resonances. However the most in- 
teresting singular power law asymptote next to the edges 
can be captured by a simpler approach which extends the 
familiar representation of electronic operators through 
Bose fields [18]. We briefly discuss it below. 

First we separate fast oscillatory modes at each edge 



out from the first band e < Epi and the second band 
Ef2 < £ < Ep3 respectively. Similarly C 2 is the sum 
over momentum of a hole particle excitations provided 
that a hole is taken out from the "gap" between Epi and 

E F2 . For example C 2 = (X)f>o ~ S^ 32 ) (coser - l)/e. 
Computing these integrals at r ^> h/\Eij\ one obtains 

Ct (r) = - log r + Si YJ Sj log \Eij I (15) 

The r-independent term in ( [T5| ) (a zero mode of the Bose 
field) explains the pre-factor in (13). Together (13) and 
dl5|) produce a canonical electronic correlation function 



In the Bose representation, Green's function ([5| is a sum 
of edge components 



c(t,x ) 



P F iX(t) 



where ipi(t) are slowly changing modes and x(t) = xo — 
vpt. Then we represent each slow modes through com- 
ponents of the Bose field as 



■£%<Pi 



(13) 



Under this representation gradients of components of the 
Bose field 



(14) 



are density of slow electronic modes. The Bose field is 
a sum of its components. The origin of the important 
factor in front of e~ £i ^ Pi will be clear in a moment. 

The component of the Bose field represent particle- 
holes excitations close to each edge. At TEij ^> H 
they can be treated as independent canonical Bose fields. 
Their variances Ci(ti,t 2 ) = -^(^| (^Pi(h) ~ Pi(h)) 2 |^) 
are not difficult to compute. As follows form ( |14| ): Ci 
are electronic density correlations at each edge:, C\ and 
C3 are sums of (coser — l)/e over all possible energies of 
particle-hole excitations provided that a particle is taken 



G{hM) = £ ^ -11 /•-•,, J e^GiihM), (17) 

Q. — ^ e (ei-a)(<^i(t2)-Vi(*i))^ J~J^ e -a(Vi(*2)-Vi(*i))^ _ 



e {\-2e % a)C %e C _ 



(18) 



We obtain Green's function in the form of two factors: 
'closed loops' e c where each edge contributes equally and 
'open lines' Li corresponding to each edge [Tj. 



Li =e i e^ 1 



C = —3a 2 log r + a 2 ^ SiSj log \E^ 



(19) 
(20) 
(21) 



This prompts the formula for the current (11). 
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